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ABSTRACT 

We present a model for launching relativistic jets in active galactic nuclei (AGN) 
from an accreting Kerr black hole (BH) as an effect of the rotation of the space-time, 
where the gravitational energy of the accretion disc inside the ergosphere, which can 
be increased by the BH rotational energy transferred to the disc via closed magnetic 
field lines that connect the BH to the disc (BH-disc magnetic connection), is converted 
into jet energy. The main role of the BH-disc magnetic connection is to provide the 
source of energy for the jets when the mass accretion rate is very low. We assume that 
the jets are launched from the disc inside the BH ergosphere, where the rotational 
effects of the space-time become much stronger, being further accelerated by magnetic 
processes. Inside the ergosphere, we consider a split topology of the magnetic field, 
where parts of the disc connect to the BH and other parts to the jets via magnetic 
field lines. The rotation of the space-time channels a fraction of the disc energy (i.e., 
the gravitational energy of the disc plus the rotational energy of the BH which is 
deposited into the disc by magnetic connection) into a population of particles that 
escape from the disc surfaces, carrying away mass, energy and angular momentum 
in the form of jets, allowing the remaining disc gas to accrete. In this picture, the 
BH can undergo recurring episodes of its activity with: (i) a first phase when the 
accretion power dominates and (ii) a second phase when the BH spin-down power 
dominates. In the limit of the spin-down power regime, the model proposed here can 
be regarded as a variant of the Blandford-Znajek mechanism, where the BH rotational 
energy is transferred to the disc inside the ergosphere and then used to drive the jets. 
As a result, the jets driven from a disc inside the BH ergosphere can have a relatively 
strong power for low mass accretion rates. We use general-relativistic conservation 
laws to calculate the mass flow rate into the jets, the launching power of the jets 
and the angular momentum transported by the jets for BHs with a spin parameter 
a* ^ 0.95. As far as the BH is concerned, it can (i) spin up by accreting matter and 
(ii) spin down due to the magnetic counter-acting torque on the BH. We found that 
a stationary state of the BH (a* = const) can be reached if the mass accretion rate 
is larger than m ^ 0.001. The maximum value of the BH spin parameter depends on 
m being less but close to 0.9982 (Thome's model). In addition, the maximum AGN 
lifetime can be much longer than ~ 10^ yr when using the BH spin-down power. 

Key words: accretion, accretion discs - black hole physics - magnetic fields - galax- 
ies :jets 



1 INTRODUCTION 

Relativistic jets are highly collimated plasma outflows 
present in extragalactic radio sources, which are associated 
with many AGN. If the matter accreted by a BH (at the cen- 
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tre of an AGN) has enough angular momentum compared 
to that of a particle moving in a circular orbit around and 
near the BH, an accretion disc can be formed. The launch- 
ing power of the AGN jet can generally be provided by the 
accretion disc, by the BH rotation, or both. Moreover, as 
the jet is launched, the BH can evolve towards a stationary 
state with a spin parameter whose maximum value is less 
but close to one (a* < 1, where —1 ^ a, ^ +!)• One can 
consider the launching power of the jet to be a fraction of 
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the disc power. A number of questions come to mind: Is this 
fraction generally valid for astrophysical jets from BHs with 
the same mass and spin? Can the disc manage to launch 
the jet by itself as the BH accretes at low rates? How does 
the magnetic field get involved? Can the BH take over and 
support the disc to launch the jet as the mass accretion rate 
goes down? How does the BH spin evolve while the jet is 
launched, and what is the maximum spin parameter in this 
case? We try to answer these questions using the model pro- 
posed in this paper. 

A super- massive BH (M ~ 10® M©) can be fed and spun 
up by accreting matter with a consistent sense of the angular 
momentu m [the first calc ulations for a Kerr BH were per- 
forme d by Bardeen (Il970|)| or by merging with another BH 



(e.g., iBerti fc Volonteril l2008l : ICergelv fc BiermannI |2009| ). 



The general relativistic effects on the structure of the inner 
regions of a n accretion disc surr o undin g a K err BH were first 
studie d by iNovikov fc Thornd (|l973h and IPage fc Thornj 
1 19741 ) using Bardeen et al.'s (1972) orthonormal frames of 
the locally non-rotating observers. These studies resulted in 
a geome trically thin, optical l y thic k accretion disc model 
[see also IShakura fc SunvaevI (| 19731 ) for a quasi-Newtonian 
approach to the description of the disc accretion on to a 
Schwarzschild BH]. This is known as the standard, thin- 
accretion disc model. The model assumes that the disc is 
quasi-Keplerian (i.e., the radial pressure is negligible and 
the radial velocity of the flow is much smaller than its az- 
imuthal component ), and extends to th e innermost stable 
(circular) orbit [see lBardeen et al] (| 19721 ') for details on the 
orbits around a Kerr BH]. The disc is driven by internal 
viscous torque, which transports the angular momentum of 
the disc outwards, allowing the disc matter to be accreted 
on to the BH. It is assumed that the torque vanishes at 
the innermost stable orbit, so that the disc matter plunges 
into the BH carrying the specific energy and angular mo- 
mentum that it has had at the innermost stable orbit. In 
the general relativistic regime, accretion on to the BH im- 
plies the conversion of the rest-mass energy of the infalling 
matter in the BH potential wall (i.e., the gravitational en- 
ergy of the accretion disc) into kinetic and thermal energy 
of the accreted mass flow. If the thermal energy is efficiently 
radiated away, the orbiting gas becomes much cooler than 
the local virial temperature, and the disc remains geomet- 
rically thin. In the inner regions of the disc, the radiation 
pressure dominates the gas pressure. The opacity is domi- 
nated by electron scattering; i.e., the photons random-walk 
before leaving the disc as they scatter off of electrons. Since 
the half-thickness of the disc at a given radius, r, is much 
smaller than the radius itself, h{r) <^ r, the disc structure 
can be described through one-dimensional (1-D) hydrody- 
nami c equations which a r e inte grated in the vertical direc- 
tion. iDonea fc BiermannI (| 19961 ) developed a model of thin 
accretion disc driven jets, which can explain the shape of 
UV spectra from an AGN when the disc is sub-Eddington. 
Models of jet/wind formation from an accretion disc typi- 
cally invoke specific magnetic field structures, as the syn- 
chrotron radio emission observed in (extra) galactic jets is 
possible only if magnetic fields are present. The jet can be 

, b y cen trifugal and 
Pavnil 19821 ). A possible 



above the disc, s hould go into unstable orbits around the BH 
(|Lvutikovll2009l ). The jets/winds can also be launched from 
either (i) a geometrically thick disc wit h, e.g., an advection- 
dominated accr etion fiow (A DAF, e.g.. iNaravan fc Yi|[l993 : 
lArmitage fc Nataraian|[l999l ). a convection-dominated flow 
(|Meier|_ I2OOII). or an advection-domi nated inflow and out- 
flow (|Blandford fc BegelmanI Il999al ), or (u) a layer lo- 
cated between the accretion disc and the BH corona 
which consists of a highly diffusive, super-Keplerian ro- 
tating and ther mally dominated by virial-hot and magne- 
tised ion-p lasma ("Huici rat et al ]|2002h . [See also the work by 
iKuncic fcB ickncU ( 200j) for the first fully analytic descrip- 
tion of a turbulent magnetohydrodynam i cs ac c retion disc 
coupled to a corona.] iFalcke fc BiermannI (| 19951 . 1 19991 ) pro- 
posed a jet-disc symbiosis model for powering jets; starting 
from the assumption that radio jets and accretion discs are 
symbiotic features present in radio quasars, these objects 
consist of a maximal jet power with a total equipartition 
(i.e., the magnetic energy fiow along the jet is comparable 
to the kinetic jet power), and the total jet power is a partic- 
ular fraction of the disc power. This fraction can be found 
by fitting the jet parameters to the observational data. 

The energy and angular momentum of a BH can be 
electromagnetically extracted in the presence of a strong 
magnetic field threading the BH and supported by ex- 
ter nal currents fiowing in th e accretion disc, as shown 
by iBlandford fc Znaielj (| 19771 ). In this case, the energy 
flux of the jets is provided by conversion of the BH ro- 
tational energy into Poynting flux, which is then dissi- 
pated a^J^xge_distaji£es_^;oin^^ BH by current instabil- 
ities (|Lvutikov fc Blandfordll2002l ). The Blandford-Znajek 
mechanism has been widely applied to jet formation in 
AGN, as well as to microquasar jets and gamma ray bursts, 
in an attempt to ma t ch a number of observational data. 
iMacdonald fc Thornd (|l982l ) explained this mechanism in 
terms of the BH membrane paradigm, in which case an imag- 
inary stretched horizon (a conducting surface located just 
outside of the BH event horizon) mimics the BH electro- 
d ynamics as seen by o utside observers. From the viewpoint 
of lThorne et~all (|l986l ) , the membrane paradigm implies not 
only the existence of a stretched horizon, but also a 3-1-1 split 
of space-time into (absolute) space and (universal) time. 
Therefore, the stretched horizon is regarded as a 2-D space- 
like surface that resides in 3-D space and ev olves in response 
to dr iving forces from the external universe (|Price fc Thornd 
Il986f ). As a result, outside observers can make measurements 
at the stretched horizon and describe the physical proper- 
ties of the horizon using pre-relativistic equations, such as 
Ohm's law. 

Different theoretical models of jet formation have 
been tested already by using numerical simulations. 
For instance, general relativistic magnetohydrodynamics 
(GRMHD) simulation results are consistent with mod- 
els of gas pressure and magnetically driven jets (e.g. , 
Koide et aT]|l999l : lMizuno et al.ll2004l:[Nishikawa et al.|[2005l : 



launched and collimate d, for instance, 



magnetic forces (e.g.. iBlandford 



condition for centrifugal launching of jets from a thin accre- 
tion disc is that the coronal particles, which are found just 



Hawlev fc Krolikll2006l : iMizuno et al.ll2007l), as weU as with 
the Blandfo r d-Znajek mechanism ( e.g., iKomissarovl [2001! : 
iKoidd I2OO3I: iMcK nnev fc Gammid |2004 ). On the other 
hand, IPunslv et al. I (|2009t ) performed fully relativistic 3-D 
MHD simulations of jets driven through the interaction of 
the magnetic field with the accreting gas in the BH ergo- 
sphere. The BH ergosphere is a part of the stationary asymp- 



Magnetic Connection Model for Launching Jets 3 



totically flat space-time (as the Kerr space-time) in which 
the Killin£| vector that corresp onds asymptotic aUy to time 
translation becomes space- like (jFriedmanlliOTSl l. Therefore, 
negative energy states of mat ter that can extract energy 
from the BH are allowed there (|Penroselll969l : lBardeen et alj 
ligyj ). The ergosphere lies outside of the BH event horizon, 
and its boundary intersects the event horizon only at the 
poles. At the stationary limit surface, an observer must move 
at the speed of light opposite the rotation of the BH just in 
order to stay still. Inside the ergosphere, the space-time itself 
is dragged in the direction of the BH rotation; i.e., nothing 
can stay there at rest with respect to distant observers, but 
it must orbit the BH in the same direction in which the BH 
rotat es. This process is c alled the dragging of inertial frames 
fe.g.. lMisner et al.lll97j ). 

The BH-disc magnetic connection, first mentione d by 
Zel'dovich & Schwartzman and quoted in lThornel(|l974l '). can 
occur and change the energy-angul ar-momentum balance 
of the accreting gas in the disc (e.g.. [ Macdonald fc Thornd 
19821: iThorrie et al.l|l98^ : lBlandfordll 19991 : Ivan Puttenll 19991 ). 



L| "(l2000al lbl. l2002al ) derived the equations for the energy 
and angular momentum transferred from a Kerr BH to a 
geometrically-thin accretion disc (which consists of a highly- 
conducting ionised gas) by magnetic connect ion, and we 
shall u se these equations. [See also the work bv lWang et al.l 
(|2003l ).] As the BH rotates relative to the disc, an electromo- 
tive force is generated. This drives a poloidal electric current 
fiowing through the BH and the disc and produces an ad- 
ditional power on the disc. From the conservation laws of 
energy and angular momentum for a thin Keplerian accre- 
tion disc torqued by a BH, (|2002ah calculated the radi- 
ation flux, the internal viscous torque and the total power 
of the disc, and f ound th at the disc can radiate even with- 
out accretion. |l| (|2002bl ) also looked for observational sig- 
natures of the BH-disc magnetic connection as more energy 
is radiated away from the disc and showed that the mag- 
netic connection can produce a very steep emiss ivity com- 

J iared to the standard, thin-accretion disc model. lUzdenskvl 
20o4 12005I ) obtained the numerical solution of the Grad- 
Shafranov equation for a BH-disc magnetic-connection con- 
figuration in the case of both Schwarzschild and Kerr BHs. 
The Grad-Shafranov equation is a non-linear, partial differ- 
ential equation that describes the magnetic fiux distribution 
of plasma in an axisymmetric system. Uzdensky found that 
this BH-disc magnetic connection can only be maintained 
very close to the BH (see in the next section). In recent 
years, a number of models that also include the BH-disc 
magnetic connection have been developed. A BH magnetic 
field configuration with both open and c losed magnetic field 
lines was considered by I Lei et al.l (|2005l '). who described the 
field configuration by the half-opening angle of the mag- 
netic flux tube on the horizon, which is determined by the 
mapping relation between the angular coordinate on the BH 
horizon and the r adial coordinated on the accretion disc. 
IWang et al.l (|2007l ') proposed a toy model for the magnetic 
connection, in which case a poloidal magnetic field is gen- 
erated by a single electric current fiowing in the equatorial 



^ Since the Kerr space-time is stationary (i.e., time-independent) 
and axially symmetric, there ar e two Killi ng vectors associated 
with these two symmetries (e.g.. [Kerrll2007l) . 



plane around a Kerr BH. lMa et al.l |20o3) derived the energy 
and angular momentum fiuxes for a Kerr BH surrounded by 
an advection-dominated accretion disc. To solve the equa- 
tions of t he accretion flow, they used a pseudo-Newtonian 
potential. iGan et al] (|2009l ) solved the dynamic equations 
for a disc-corona system and sim ulated its X-ray s pectra 
by using the Monte Carlo method. IZhao et al.l (|2009l ) stud- 
ied the magnetic field configuration generated by a toroidal 
distributed continuously in a thin accretion disc, as well as 
the role of magnetic reconnection in the disc to produce 
quasi-peri odic oscillati o ns in BH binaries. In the context of 
GRMHD. lKoide et all (|2006l ) presented a 2-D GRMHD re- 
sult of jet formation driven by a magnetic field produced 
by a current loop near a rapidly-rotating BH, in which case 
the magnetic fiux tubes connect the region between the BH 
ergosphere and a co-rotating accretion disc. Furthermore, 
relativistic Poynting jets driven from the inner region of an 
accretion disc that is initiall y threaded by a dipole-li k e mag - 
netic field were studied by iLovelace fc Romanoval (|2003h . 
Their model is derived from the special relativi stic equation 
for a force-free electromagnetic field. [See also iLvnden-Bell 
(1996, 200I).] 

In this paper, we propose a model for launching rela- 
tivistic jets from a (geometrically-thin) disc inside the er- 
gosphere as an effect of the rotation of the space-time. We 
consider here the BH-disc magnetic connection, whose main 
role is to provide the source of energy for the jets when 
the mass accretion rate is very low. We use the general rel- 
ativistic form of the conservation laws for the matter in 
a thin accretion disc to describe the disc structure when 
both the BH-disc magnetic connection and the jet forma- 
ti on are considered. The model is based on the calc ulations 
oflNovikov fc Thornd (| 19731 ). I Page fc Thornd (| 1 9741 ) and Q 
('2002aD, be ing mainly influenced by th e work of IZnaiekl 
tl97i ) and iMacdonald fc Thornd (| 19821 ). [Some incipient 
ideas which are at t he base of this model were exposed in 
iDutan fc BiermannI (|2005al lbl).] This is the first work that 
studies the process of jet launching from a geometrically- 
thin accretion disc inside the BH ergosphere when the en- 
ergy and angular momentum are transferred from the BH 
to this region of the accretion disc via closed magnetic field 
lines, within the framework of general relativity. An impor- 
tant result of the model, with impact on observation of the 
AGN jets, is that the power of the jets does not depend 
linearly on the mass accretion rate all the way down to 
very low accretion rates for BHs of a given m ass. This re- 
sult is different from that of lAUen et al.l (|2006l ). who found 
a linear dependence between the power of the jet and the 
mass accretion rate by considering a spherical Bondi-type 
accretion on to BHs (in which case the accreting matter 
has zero or very low angular momentum). In their calcu- 
lations, the power of the jet is estimated from the energy 
and time scale required to inflate the cavity observed in 
the surrounding X-ray emitting gas. The model proposed 
here combines two regimes associated with the driving of 
the jets, an accretion power regime and a (BH) spin-down 
power regime, where the switch from the former to the latter 
regime corresponds to a mass accretion rate of m ~ 10~^ *. 
In the accretion power regime, the power of the jets is lin- 
early dependent on the mass accretion rate, whereas in the 
spin-down power regime the power of the jets depends very 
weakly on the mass accretion rate. In the accretion power 
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regime, the energy and angular momentum are extracted 
and transported away from the disc inside the BH ergo- 
sphere by both the kinetic fiux of particles and Poynting 
flux in the form of jets. Instead, in the spin-down power 
regime the energy and angular momentum are extracted and 
carried away from the disc inside the ergosphere predomi- 
nantly in the form of Poynting flux with just little amount 
of kinetic flux of particles . The work p r esent ed in this paper 
is different from that of IWang et al.l (|2008l ) , in which the 
production of Poynting flux jets is associated with a com- 
bination of the Blandford-Znajek mechanism, the BH-disc 
connection and the Blandford-Payne mechanism. Further- 
more, we argue that the accretion, which is initially at either 
close to the Eddington rate or at low rates, can be driven 
in a non- radiant, geometrically thin and quasi-Keplerian 
disc inside the BH ergosphere by the external jet torque. 
This is distinctly different from opti cally thin, advection - 
dominated accretion flow models ('e.g.. lNaravan fc Yilll994 ). 
in which the accretion at low rates is advection-dominated 
(i.e., the thermal energy generated via viscous dissipation is 
mostly retained by the accreted mass flow rather than be- 
ing radiated, and the energy is advected in towards the BH) 
and quasi-radial (i.e., the flow is geometrically thick, having 
roughly a spherical structure). Our first goal in the present 
work is to obtain estimations for jet-related quantities, in 
particular the mass loading, power and Lorentz factor of 
the jet, when the BH rotational energy is transferred to the 
disc inside the ergosphere, and then to compare them with 
those derived from the Blandford-Znajek mechanism. In the 
limit of the spin-down power regime, the model proposed 
here can be regarded as a variant of the Blandford-Znajek 
mechanism, where the BH rotational energy is transferred 
to the disc inside the ergosphere rather than transported to 
remote astrophysical loads. Our second goal is to determine 
the upper limit of the spin parameter attained by a station- 
ary Kerr BH when both jet formation and BH-disc magnetic 
connection are considered and to investigate how the value 
of the mass accretion rate can influence the departure of the 
BH spin parameter from its theoretical maximum limit of 
a, = 1. 

In Section (2] we describe the assumptions of the model. 
In Section O we derive the mass flow rate into the jets. 
Using the general-relativistic conservation laws for matter 
in the accretion disc (Section [4|, we derive the launching 
power of the jets (Section [S]) and the angular momentum 
removed by the jets (Section [6]). In Section [71 we calculate 
the efficiency of launching the jets and show that when the 
BH accretes at low rates, the spin-down of the BH is an 
efficient mechanism to launch the jets via the accretion disc. 
In Section |8l we study the spin evolution of the BH and 
discuss conditions of BH stationary states for given mass 
accretion rates. In Section [9] we refer to the long lifetime 
of an AGN from the BH spin-down power as a particular 
relevance of the proposed model to the observational data. 
In Section If 01 we present a summary of the key points, as 
well as our conclusions, and suggest further work to be done. 



2 BASIC ASSUMPTIONS 

• We consider that matter outside of the BH has negli- 
gible gravitational effects compared to the BH gravity and 




spin parameteia« 



Figure 1. Schematic representation of the inner part of the ac- 
cretion disc-BH-je t system, where the BH is represented by the 
stretched horizon jThorne et al.lll986l) . Above the surface of the 
disc inside the BH ergosphere, the closed magnetic field lines 
(solid lines) do not cross the open magnetic field lines (dashed 
lines); they overlap only in the line-of-sight projection. For some 
explanation on the structure of the magnetic field in the disc in- 
side the BH ergosphere, the reader is referred to the text below. 



that an accretion disc settles down in the equatorial plane 
of a Kerr BH. We assume that the accretion disc is ge- 
ometrically thin and quasi-Keplerian, therefore the physi- 
cal quantities integrated over the vertical direction can be 
used in order to describe the radial structure of the disc. 
The inner part of the accretion disc is located inside the 
BH ergosphere, extending from the stationary limit surface 
(or ergosphere) inward to the innermost stable orbit (see 
Fig. [T]). In the BH equatorial plane, the stationary limit 
surface is located at r^i = 2rg and does not depend on 
the BH spin parameter, whereas the radius of t he inner- 
most stable orbit depends on o* [cf. eq. 2.21 of iBardeenI 
(,f970. )]: so that, rmsia.*) = 2rg for a* = 0.95 and ~ f-2rg 
for a* ~ f. Here, rg = GM/c? is the gravitational radius 
[rg = rJ(M/fO^MQ) = f.48 x f 0'*(M/f O^Mq) cm], G is 
the Newtonian gravitational constant, M is the BH mass 
and c is the speed of light. 

• We adopt th e viewp oint of the BH membrane paradigm 
of [rh orne et ahl l)f 986h and represent the BH by the 
stretched horizon. 

• We consider the case of rapidly-spinning BHs with 
a spin parame ter a* ^ 0.95, based on the argument by 
iBardeenI l) 19731 ) that a strong preference for a particle to 
orbit in the equatorial plane requires the BH spin parame- 
ter to be close to its maximum value. 

• We consider that closed ma gnetic field lines con nect th e 
BH to the accretion disc (e.g.. iBlandfordllf 9991 : [Lil l2000bl ). 
The poloidal component of the magnetic field transfers an- 
gular momentum and energy (in the form of Poynting flux) 
from the BH to the disc, thereby increasing the amount of 
the gravitational energy which is released from the accretion 
disc. This energy is liberated very close to the BH, where 
most of the gravitational energy of the accretion disc is avail- 
able (in our case, from the disc inside the BH ergosphere). 
Next, we assume that the disc energy (which, once again. 
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represents the disc gravitational energy plus the BH rota- 
tional energy deposited into the disc by magnetic connec- 
tion) is used to launch the jets. In this way, the disc remains 
cool and geometrically thin. By comparison, in a standard, 
thin accretion disc, the thermal energy is efficiently radiated 
away, so that the orbiting gas becomes much cooler than the 
local virial temperature, and the disc remains geometrically 
thin. 

• As the numerical results obtained bv lUzdenskvl (|2005l ') 
indicate, a closed magnetic field configuration can only be 
maintained in a region close to the BH. This region is limited 
to the radius of ~ 12 rms = 72 rg for a static BH (a, = 0) 
and decreases to ^ 6 rms ~ 20.4 rg as the BH spin param- 
eter increases to a, — 0.7. Therefore, it might be possible 
that, for a BH with a, ^ 0.95, the closed magnetic field 
configuration to exist only in the region bounded by the BH 
ergosphere. 

• We constrain the jet formation to the disc inside the 
ergosphere, where the rotational effects of the space-time 
become much stronger. The slope of the specific energy of 
disc particles (Eq. I13|) for a, ^ 0.95 steepens in the BH er- 
gosphere, indicating that most of the gravitational energy 
of the accretion disc is released from there. Moreover, inside 
the ergosphere, the space-time is dragged in the direction of 
the BH rotation and everything must co-rotate with the BH. 
(We shall ignore the possibility for the plasma inside the BH 
ergosphere to have negative mechanical energy at infinity.) 
Outside the ergosphere, it is plausible that patches of the 
disc gas may counter-rotate with the BH due to some insta- 
bility or turbulence in the disc. If some magnetic field lines 
are frozen in such patches, they likely reconnect with the co- 
rotating magnetic field lines, and the closed magnetic field 
configuration can be destroyed. On the other hand, following 
the discussion above, it would appear that a closed magnetic 
field configuration will exist only in the region bounded by 
the BH ergosphere. And it is this magnetic field that trans- 
fers rotational energy from the BH to the accretion disc. 

• A more general magnetic field threading the BH would 
consist of a combination of closed magnetic field lines, as 
discussed above, and open magnetic field lines (extending 
from the BH to infinity) , which can enable the extraction of 
BH rotational energy through the Blandford-Znajek mech- 
anism. Now, one can wonder: (i) How much of the rotational 
energy of the BH is directly extracted by jets and how much 
is transferred to the disc, and (ii) under which specific con- 
ditions? These questions can be addressed in some future 
work. We limit the present work to the case of a BH threaded 
only by closed magnetic field lines and compare the values 
of some basic physical quantities (such as the power) of a 
jet driven solely from the disc inside the BH ergosphere with 
those of a jet powered by the Blandford-Znajek mechanism. 

• We assume that the jets remove angular momentum 
from the disc, enabling the accretion process. The external 
torques acting on the disc inside the BH ergosphere (i.e., 
magnetic BH-disc torque and jet torque) can dominate ove r 
the internal viscous torque of the disc (|Blandfordll200lf l. 
where the internal torque can be produced by m agneto- 
rotational instability (e.g.. [Balbus fc HawlevI Il99ll ). which 
is driven by the free energy available from the differential 
rotation of the gas fiow. 

• To form the jet, at some point, the particles must cross 
the magnetic field lines. One can picture this as being due 



to dr i fts and inst abilities inside the disc (|Balbus fc HawlevI 
1 1994 Il996l . Il998h . In addition to the closed magnetic field 
lines that connect the BH to the disc inside the ergosphere, 
the magnetic field can be represented similar to the mag- 
netic field on the Sun surface. In a direct top-down view, it 
looks mottled, with arcs of magnetic fields connecting dif- 
ferent regions and other magnetic flux tubes extending into 
fr ee space and allowin g matter to flow out [see also fig. 36 
of iThorne et al.l (|l986l l]. [An accretion disc with a magnetic 
fie ld configura t ion si milar to that of the Sun was studied 
bv lYuan etal] (|2009l ') for driving episodic magnetohydrody- 
namic (MHD) jets.] However, the closed (BH-disc) and open 
(disc-infinity) magnetic field lines are supposed to be domi- 
nant in the model proposed here. 

• To examine the fate of the magn etic field in the space 
around the BH, we follow the ideas of iThoriie et all (|1986[ ) 
on the 'cleaning' of the magnetic field by the BH stretched 
horizon to maintain an ordered magnetic field [see fig. 36 of 
iThorne et al] (|l986l )]. Since the magnetic field is frozen into 
the accretion disc, the field lines are transported in toward 
the BH by the accretion flow. Once the flow reaches the 
innermost stable orbit, it drops out of the disc and falls di- 
rectly into the BH, and becomes causally disconnected from 
the field lines to which it was attached. However, the fiux 
conservation law assures that the field lines, although dis- 
connected from their source, will be pushed onto the BH by 
the Maxwell pressure of the adjacent field lines or thread 
through the near-vacuum region between the accretion disc 
and the BH or pushed back into the disc via Rayleigh- Taylor 
instabilities. If this cleaning works efficiently enough, an ap- 
proximately stationary and axisymmetric magnetic field of 
the type de scribed in the paragraph above can be formed 
see alsoO l|2000j )]. otherwise the process discussed in this 
paper is not continuous and a different description than 
that presented here must be employed. Furthermore, our ap- 
proach here is simple in the sense that with the assumption 
of an approximately stationary and axisymmetric magnetic 
field, the outflow near the launching region is thereby taken 
to be quasi-stationary {d/dt = 0). Therefore, we are lim- 
ited by not being able to specify time-dependent effects in 
the outflow. We concentrate here on the launching region of 
the jets, determining the mass flow rate into the jets and 
other conserved quantities of the outflow, and discuss the 
importance of the magnetic field in accelerating the jets on 
sub-parsec scales without solving the entire fiow problem 
(see Section [5]). 

• To allow for jet formation, we assume that the surface 
density of the disc and the mean radial velocity combine 
themselves in such a way to keep their product constant 
over the disc inside the BH ergosphere (see Section |3]). This 
can be thought of as being an effect of the BH ro t ation . 
By comparison, in the model of iBlandford fc Pavnel (|l982l ) 
the disc particles are driven upwards by the gradient of the 
pressure in the disc to fill the corona around the disc. The 
plasma should be sufficient to produce whatever charge and 
current densities are required for an MHD flow. The mag- 
netic fleld lines that pass through the corona will be bent 
from near vertical to make a certain angle with the disc 
surface. If this angle is greater than 60°, the particles will 
fling outwards under centrifugal forces. Then, as magnetic 
stresses become important, the particles will be further ac- 
celerated by the gradient of the magnetic pressure. Far from 
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the disc, the particles inertia will cause the magnetic field 
to become increasingly toroidal. As a result, the flow will 
be collimated by the magnetic hoop stresses to a cylindrical 
shape outside the (outer) light cylinder. In the model pro- 
posed here, the rotation of the space-time is responsible for 
ejecting particles in the direction perpendicular to the disc. 
The escape particles then slide along the open magnetic field 
lines (dashed lines in Fig. [1]) , being accelerated by magnetic 
forces (see Section [5]). 

• Due to the motion of the disc plasma in which mag- 
netic field lines are frozen, it is also possible that magnetic 
reconnection will take place at the interface of closed and 
open magnetic field lines. (Such a non-ideal MHD effect, 
which is common to impulsive and gradual fiares, can con- 
vert magnetic energy into kinetic energy and thermal en- 
ergy, and topologically change field lines; this can have sig- 
nificant consequences for the global evolution of a system. 
This effect is also connected with a violation of the magnetic 
ffux conservation.) Such a magnetic reconnection that allows 
outflows to be driven has been observed in time-dependent 
resistive MHD simulatio ns of jet /wind formation from neu- 
tron stars performed by iRomanova et all Hooi). This role 
of magnetic reconnection in driving outflows can also be ex- 
tended to accreting BH systems. To estimate the rate of the 
magnetic energy that can be extracted through reconnec- 
tion, we follow the work by de Gouveia dal Pino fc LazarianI 
l|2005l ) ; lde"Gouveia Dal Pino et all (l2010l'). If we applv eo. 
12 of IdT Gouveia dal Pino fc LazarianI (|2005h to a BH of 
10^ Mq, we obtain a value of ~ 10*'^ erg s~^, which is about 
two orders of magnitude smaller than the minimum value of 
the power of the jet that we obtain in the model proposed 
here. 

• We use the Kerr metric (|Kerilll963l ) in cylindrical coor- 
dinates. In and near the BH equatorial plane, the metric is 
given by 



ds 



-e at + e 



2 ,p , 



- bjdt) + e ^dr -I- dz 



(1) 



where r, (f) and z are defined as the cylindrical coordinates 
in the asymptotic rest frame, a nd t is the physical ti me of 
an observer removed to infinity (jPaee fc Thornelll974h . The 
components of the metric tensor in ^ are specified by 



A ' 



2fi 



where the metric functions read 



2rgr + , A = r"^ + r^a' -f Irf/ra' 



(2) 



(3) 



where a = J/{Mc) is the angular momentum of the BH 
about the spinning axis (J) per unit mass and per speed of 
light. The BH spin parameter is defined as a, = J/Jmax [= 
a/vg), where Jmax = GM^/c is the maximal angular mo- 
mentum of the BH. 



3 MASS FLOW RATE INTO THE JETS 

First, let us explain the terms used for the (rest) mass rates 
of the flow, which are measured by observers at inflnity. 
The term mass accretion rate refers to the mass flow rate 
through the disc up to the ergosphere (Md = rhMEdd)- In 
the ergosphere. Ado is divided into the mass outflow rate 
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'jet 



'ms 



innermost 
stable orbit 



Msl 



stationary 
limit surface 



Figure 2. Schematic representation of the mass flow through the 
disc inside the BH ergosphere. 



into the jets (A/jcts) and the mass inflow rate on to the BH 
(Min). Now, we can write the mass outflow rate as 



Mjots = Min(rsl) - Min(rms), 



(4) 



where Min(rsi) and Min(rinB) denote the mass inflow rate 
at the stationary limit surface and at the innermost stable 
orbit, respectively. The mass outflow rate can also be ex- 
pressed as 



(5) 



where the parameter gjots indicates the fraction of the mass 
of the disc inside the BH ergosphere that goes into the jets. 
Figure [5] shows a schematic representation of the mass flow 
in the disc inside the ergosphere. 

The amount of mass that flows inward across a cylin- 
der of radius r during a co ordinate time interval At , when 
averaged by the method in lNovikov fc Thornd (|l973l ). is 



M : 



1/2 



(6) 



where y^'Jgl = e"^'''"'''^ = r is the square root of the metric 
determinant (Eqs. [T] and (2)1, E = 2/i < po > is the surface 
density of the disc (with h being the half-thickness of the disc 
and < po > the density of rest mass), is the mass-averaged 
radial velocity and T) — (1 — 2/r»-|- a1/r1) is one of the 
functions used to calculate general relativistic corrections to 
the Newtonian accretion disc structure. Here, r, = r/r-g is 
the dimensionless radius. 

Next, we estimate the mass inflow rate at one speciflc 
radius of the disc inside the BH ergosphere by using so 
that the mass outflow rate becomes 



Mjets = [-27rr- E i;^ - [-27rr- E i;" . (7) 

Now, the mass inflow rate at the stationary limit surface 
is also given by mAfEdd, and then 



Afjets = mAfedd 



Af(rms) 



A^(rei) 



(8) 



By analogy w ith a standard, thin accre tion disc [eqs. 
5.9.5 and 5.9.10 of iNovikov fc Thornd (|l973l ')]. the product 
of the surface density and the mass-averaged radial veloc- 
ity can depend on the disc radius as Eu"" oc . One can 
recover a standard, thin accretion disc by setting p = — 1, 
in which case the ma s s accr etion rate [Eq. [6l also eq. 5.6.2 
of iNovikov fc Thornd (|l973l )] is independent of the radius. 
Now, instead of trying to obtain an exact condition for jet 
launching, let us simply take p = 0, as this gets rid of the 
requirement to know E and v"^ precisely. Let us examine 
the consequences of this choice. Clearly, p = corresponds 
to a product Ew"" constant for any radius of the disc in- 
side the ergosphere, and consequently the mass inflow rate 
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Figure 3. Mass outflow parameter (gjcts) ^ a- function of the BH 
spin parameter (a*). For a* ~ 1, almost whole material of the disc 
inside the BH crgospherc flows into the jets (gjcts — 0.98); that 
is, the BH almost stops being fed by accreting matter. 



on to the BH decreases with radius as Min(r) ~ r, and 
the difference flows into the jets. Then, Min(r) at a given 
radius is specified by the general relativistic factor D^^^, 
which depends on the BH spin parameter. For an ADAF 
disc that allows outfiows, Mi n (r) ^ r '^, where ^ q < 1 [eq. 
13 o f iBlandford fc Beeelmanl (|l999bl )]. Therefore, our choice 
of q = 1 (or p = ) repres e nts a l imiting case in the model 
bv lBlandford fc Begelnianl (| 1999^ 1 . 

Using equations ((5]) and ([8]) , as well as the expression of 
we obtain the fraction of the mass inflow that goes into 
the jets (or the mass outfiow parameter) as 



(Jjcts(ffl*) 



1/2 



(9) 



where rms» = rms/rg and r^i, = r^i/rg. 

FigureOshows the mass outflow parameter as a function 
of the BH spin parameter. For a* — 0.95, the mass outflow 
into the jets is only about 8 per cent of the available mass 
inflow through the disc inside the ergosphere. Instead, for 
a spin parameter near the maximal one (a* ~ 1), the mass 
outflow increases to about 98 per cent of the available mass 
inflow. This means that in the case of near maximal spin, 
the BH almost stops being fed by accreting matter. 

Suppose the inner disc would have been extended be- 
yond the stationary limit surface. In this case, the disc par- 
ticles can form the jets (gjots > 0) if and only if the BH spin 
parameter were a» > 0.755 (plot not shown). 

We mention that the results presented in this section 
are valid for our choice of p = 0. This, of course, need not 
be a necessary condition for jet launching, since we have ex- 
amined the mass outflow parameter for just one value of p, 
that which makes T.v^ constant for any radius of the disc 
inside the BH ergosphere and gjets dependent only on the 
BH spin parameter, but it is certainly sufficient. The BH 
rotation causes an outflow of particles from the disc, where 
the energy (and angular momentum) carried by the escape 
particles is taken from the accretion disc. The escape par- 
ticles then slide along the open magnetic field lines, being 
accelerated by magnetic forces (see Section [5]). 



4 ANGULAR MOMENTUM AND ENERGY 
CONSERVATION LAWS 

To describe the structure of the disc inside the BH ergo- 
sphere, we use the ang ular momentum and e nergy conser- 
vation laws derived by IPage fc Thornj (|l974h and include 
both the BH-disc magnetic connection an d the jet forma- 
tion. W hen deriving the conservation laws, IPage fc Thornd 
do not make any assumption about the type of stress- 
energy present (e.g., magnetic fi elds, viscous stresses, e tc.). 
[The calculations performed by IPage fc Thornd (| 19741 ) are 
valid even if the disc is highly dynamical, but can also be ap- 
plied to steady-state and quasi-steady-state discs, in which 
case the mass accretion rate is constant throughout the disc] 
Here, we consider that the removal of the angular momen- 
tum of the disc inside the ergosphere can be produced by the 
external jet torque and that the external torques acting on 
the disc inside the ergosphere (i.e., BH-disc magnetic torque 
and jet to rque) dominate over the internal viscous torque of 
the disc (|Blandfordll200ll ). In this case, we can write the 
angular momentum conservation law0 as 



d_ 



1 - gjcts) MocL' + AixrH = AurJV 



(10) 



where on the left-hand side, the first term describes the an- 
gular momentum carried by the accreting mass of the disc 
inside the BH ergosphere, and the second term describes 
the angular momentum transferred from the BH to the disc 
inside the ergosphere. The term on the right describes the 
angular momentum carried away by the jets. is the spe- 
cific angular momentum of a gas particle orbiting in the 
accretion disc, J is the total flux of energy (of particle and 
electromagnetic origin) carried away by jets and H is the 
flux of angular momentum transferred from the BH to the 
disc inside the ergosphere. H is defined through the mag- 
netic torque produced by th e BH on both surfaces of the 
accretion disc Thd (|Lil l2002al ) 



2-KrHdr, 



(11) 



where the limits of integration are two radii of the accretion 
disc with ri < r2. 

Similar to the angular momentum conservation law, we 
can write the energy conservation law as 



5 I" 



Mdc^BM + ^TTrHno = 4TvrJE\ (12) 



where on the left hand-side, the first term describes the rate 
of the energy fiow through the disc inside the BH ergosphere, 
and the second term is the rate at which the magnetic torque 
per unit area of the disc does work, Thd^Id (here S7d is 
the Keplerian angular velocity of the gas particles in the 
disc). The third term describes the energy flow along the 
jets. is the specific energy of a gas particle having mass 
^l and orbiting in th e same direction as the BH rotation 
(|Bardeen et aLill972h : 



E 



2rgr^''^ + rV^a 



1/2 ■ 



(13) 



The equation describing the angular momentum conservation 
is derived in Appendix lAl 



8 loana Du^an 



The flux of angular momentum transferred from the BH 
to the disc inside the ergo sphere by magnetic connection has 
the following expression (P l2002al ): 



H 



STT^r \cdr J {—dRu/dr)^ 



(14) 



where ^Jd is the flux of the poloidal magnetic field lines 
which thread the surface of the disc inside the BH ergosphere 
and fin is the BH angular velocity. The derivation of Eq. H14[) 
is based on the supposition that the accretion disc consists 
of a highly conducting ionised gas. This implies that (i) the 
accretion disc resistance is neglected in comparison with the 
BH surface resistance and (ii) the magnetic field lines are 
frozen in the accretion disc, being transported by the disc 
gas and rotating with fio. On the other hand, the angular 
velocity of the magnetic field lines threading the horizon 
is JIh, due to the effect of the frame-dragging at the BH 
horizon. For a* > 0.35 and r ^ rms, S^h > ^o, so that the 
BH transfers energy (and angular momentum) to the disc. 
For a, < 0.35, flu < fio, and this time the accretion disc 
transfers energy (and angular momentum) to the BH. For 
a* = 0.35, Qh = flo; this condition implies that there is no 
energy (nor angular momentum) transfer between the BH 
and the accretion disc by magnetic connection. 



5 LAUNCHING POWER OF THE JETS 

We are now in the position to calculate the launching power 
of the jets with the help of the conservation laws previously 
derived. First, we define the launching power of both jets as 



-jets 



2 / " 2TiJEU-dr. 



(15) 



Integrating the equation of the energy conservation law (Eq. 
I12p over the disc inside the BH ergosphere, we find the 
launching power of the jets. 



- jets 



(1 - fto 



I Muc ( El - El, 



A-K I rHQ,udr. 



(16) 

The first term describes the rest energy of the accreting 
matter on to the BH and the second term describes the 
energy transfer from the rotating BH to the disc inside the 
ergosphere. and are the specific energy of the gas 
particle (Eq. I13p evaluated at the stationary limit surface 
and at the innermost stable orbit, respectively. 

Using Eq. (|14|l . we obtain the launching power of the 
jets as 



Pjots = (1 - ftets) MUC {eI - El 

Oh - 



+ 



1 /dl-D 

2^ 



Or 



cdr J {—dRn/dr) 



QTjdr, 



(17) 



where the angular velocities of the BH and the accretion 
disc, respectively, are 



2^g l + (l-a2)i/2 



(18) 



to evaluate both ^'d and {—dRu/dr). First, we write the 
magnetic fiux that threads the accretion disc surface. 



Bu{dS), 



(20) 



where Bb is the poloidal component of the magnetic field 
that threads the disc. The surface area between two equa- 
torial surfaces in a Kerr space-time can be calculated from 



(d5')2=o = det dr d(j), 
where the determinant of the surface metric is 
det 3(^0) = 





gr<t> 




g2M 





g<pr 


94"t> 







e2V, 



A 
A' 



(21) 



(22) 



This result follows from Eqs. d), ([2]), and With these, 
the surface area in Eq. (|2ip reads 



(d5).=o = ( ^ 



1/2 



27r dr. 



(23) 



The poloidal component of the magnetic field that 
threads the BH horizon Bu and the poloidal component of 
the magnetic field at the inner edge of the accretion disc 
Bu{rms) can be of the same order (e.g., iLivio et aLlll999( ) 
and related by 



-Bh = CBoirms) , where C > 1. 



(24) 



On the other hand, the poloidal component of the magnetic 
field that threads the accretion disc su rface scales as Bu oc 
r~", where < n < 3 (|Blandfordlll976l 'l. Consequently, 



Bu = Buirms) 



Bh 

c 



(25) 



Since the BH horizon behav es, in some aspects, like a ro- 
tating condu c ting surface (e.g.. [Carted Il973l : iDamourl Il978l : 
IZnaiek|[l97i : iThorne et all 1 19861 '). itcan be thought of as 
being a 'battery' driving currents around a c ircuit. The en - 
ergy for this comes from the BH rotation (jZnaiekl Il978l l. 
The internal resistance of the battery in the horizon, i.e., 
the resistance between two magnetic surfaces that thread 
the horizon, is 



dRu — Rh 



(26) 



where Ru = 4n/c = 377 ohm, dl is the horizon distance be- 
tween two magnetic surfaces (see Fig. [T]), 2nrH is the cylin- 
drical circumference at r = rn and rn = rg[l-|-(l- 



2n1/2i 



Tgru* is the radius of the BH horizon (|Thorne et al.lll986r i. 

The voltage difference generated by the BH has a max- 
imum magnitude of V = nn^&H, where ^'h = BuAh is the 
magnetic fiux threading the BH and An = 87rrgrH is the 
surface area of the BH. Assuming that the magnetic field is 
carried into the BH by the accreted disc gas, we set the BH 
potential drop to the energy of the gas particles carried into 
the BH, the latter being the particle specific energy at the 
innermost stable orbit. Suppose that during a first epoch, 
the BH accretes at a rate approximately equal to the Ed- 
dington rate|f] This supposition provides — M^ccElsC^ . 



fir 



3/2 



■ flu*- 



(19) 



To calculate the launching power of the jets, we need 



^ The Eddington accretion rate is defined from the Eddington 
luminosity as M^dd = ^Edd/i^c^) = 47rGM/(eKxc), where e is 
the efficiency of converting the gravitational energy of the accre- 
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Therefore, the maximum value of the magnetic field that 
threads the BH horizon is 



(27) 



where a^jim = 0.9982 is the BH limiti ng spin in th e case 
of a radiatively efficient accretion disc (|Thornelll973 ). and 
the corresponding particle specific energy at the innermost 
stable orbit is = 0.6759. Although this limit of the 

BH spin may be even closer to its maximum value a* ~ f , 
it produces a negligible variation in the maximum value of 
the BH magnetic field. Using the expression of the gravita- 
tional radius, the maximum value of the magnetic field that 
threads the BH horizon fEg. I27() becomes 



Ar rE' 



M 



Bh"" = 0.56 X fO-* 



/ M 



-1/2 



gauss. 



(28) 



(29) 



This r esult is simi lar to the c alcula tion performed bv lZnaie^ 
(jl978h . [See also iLovelacd (jl976h .] The only difference is 
that we set the BH potential drop to the specific e nergy of 
the pa rticles at the innermost stable orbit, whereas IZnaiekl 
l|l978h makes use of the fact that the Eddington luminosity 
sets an upper bound on the radiation pressure (as the disc 
is radiatively efficient), and thus V'^ ~ -J/Edd. 

The continuum of the magnetic field within a narrow 
strip between two magnetic surfaces, which connect the 
BH to the disc i nside the ergosphere, ci^'H = d*D (e.g., 
IWang et al.ll2007h . gives 



BH27rrH dl = -Bd ( ^ 



1/2 



2-K dr. 



Making use of Eqs. (f25ll . pejl and pO)) . we obtain 



i-dRn/dr) = ^ ^ 



2 I I r 



A 



1/2 



(30) 



(31) 



Next, we write the particle specific energy fEq. If 3|l us- 
ing the dimensionless radius and the spin parameter. Then, 
we substitute this equation (evaluated at Vsi, and rms* , re- 
spectively) together with ^ ~ ((27]), (O and ((31)) for Eq. 
(|f 7p . So, the launching power of the jets becomes 



Dmax 



M 



(lO9Af0 



(32) 



tion disc into radiation and K-p denotes the Tho mson opacity, e 
depends on the BH spin parameter as e = 1 — Ems l lThornelll974) . 
so that e = 0.06 for a Schwarzschild BH and e = 0.42 for an ex- 
tremely spinning Kerr BH. We scale the BH mass to IO^Mq, so 



that A^Edd = M],..e-^[MI\Q'^MQ), where M\ 



1.38 X lO"' 




Figure 4. Launching power of the jets as a function of the mass 
accretion rate m (Eq. I32|l for a given BH spin parameter, a* = 
0.99. The switch from an accretion power regime to a spin-down 
power regime corresponds to a mass accretion rate of m ~ 10~^'^. 



where 



C, = 



47r(^(a, 



lim)^eiin 



R. = 



1 + ajrj^ -f 2a>r^ 
1 _ 2^-1 -f air.-' ' 



1 o -1 1 -3/2 

1 — 2r, + a.r* ' 



3rr^ + 2a*r, 



3/2 



1/2 



(33) 



(34) 



For the following calculations, we consider the strength of 
the magnetic field in Eq. (|32p to be as high as its maximum 
value Bh ~ B^^^ . On the right-hand side, the first term 
represents the accretion power of the disc inside the BH er- 
gosphere and the second term represents the BH spin-down 
power transferred to the disc by magnetic connection. So, 
Eq. (|32|) can also read: 



p _ pace . pre 
jets — -'jets ~r J jc 



(35) 



g s 



In Fig. 21 we plot the launching power of the jets as a 
function of the mass accretion rate. The plot shows that a 
transition from an accretion power regime to a spin-down 
power regime is produced for m ~ fO"^ *. So, we have: (1) 
an accretion power regime in which case m > f O"^ * and the 
dominating term in the launching power of the jets is Pjcts) 
and (2) a spin-down power regime in which case m < 10~^'* 
and the dominating term in the launching power of the jets 

-Pjets- 

In Eq. H32p . the launching power of the jets depends on: 
(i) the mass accretion rate rh, (ii) the BH mass M, (iii) the 
BH spin parameter a*, (iv) the power-law index n, and (v) 
the ratio of the magnetic field strengths i^. We chose the last 
two parameters as follows: the power-law index n i s taken 
to be '2' as for a frozen magnetic field (Alfvcn f 9631 ). and C, 
is set by taking its value corresponding to the maximum of 
the launching power of the jet, which is one. Therefore, for 
the following calculations, we consider n = 2 and(" = 1. 

In Fig. [SI we plot the launching power of the jets 
as a function of the BH spin parameter, for a BH mass 
of 1O'A/0, given four values of the mass accretion rate 
(m = 1, 0.5, 0.1 and 0.01), as well as the BH spin-down power 
contribution to the jets power (bottom curve). Since the area 
of the disc inside the ergosphere increases with an increase 
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0.99 
spin parameter a. 

Figure 5. Launching power of the jets as a function of a» (Ea. l32ll 
for four values of the mass accretion rate. The mass accretion rates 
from top to bottom are 1, 0.5, O.f and O.Of. The bottom curve 
represents the power of the jets given by the BH spin-down ^j™!- 
Note that there is a sUght difference between the two bottom 
curves (red and turquoise curves). In the case of very low mass 
accretion rates, rh < O.Of, Pjcts is approximately equal to the BH 
spin-down power. 
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Figure 6. Initial magnetization parameter of the jet plasma as a 
function of the disc radius for three values of the mass accretion 
rate when the BH spin parameter is a* = 0.99. The mass accretion 
rates from top to bottom are O.OOi, O.Of, O.f and 1. For cr ^ 1, the 
jets are Poynting fiux-dominatcd outflows; i.e., the energy content 
of the jets is mainly in the magnetic part. For cr < 1, the jets are 
particle-dominated outflows. Note that ctq ~ m^^. 



of at, there is a dominating trend of the jet power to in- 
crease as well, except for a« close to the maximal value and 
rh > 0.1 where the turn-over of the curve is produced due 
to the general relativistic factor that appears in the term 
(1 — gjets) of the accretion power. In the case of the spin- 
down power regime, the jet power is ~ 10*^ erg s~^, which 
is only 10~^ of the Eddington luminosity of a 10^ solar mass 
BH. This value of Pj°ts is comparable to the maximum rate 
of energy extraction by the Blandford-Znajek mechanism, 
which is ~ 10*^ erg s~^ for a BH mass of 10^ Mp) and a* close 
to the maximal value [Eq. 4.50 of lThorne et al.l (j 19861 )]. For 
a lower mass of the BH, the jet power decreases, as the 
launching power of the jets is proportional to the BH mass. 

On sub-parsec scales, the jets are likely to be domi- 
nated by electromagnetic processes (MHD or pure electro- 
dynamic), where the energy is transported along the jets 
via Poynting flux, and are pote ntially unsta ble if significant 
thermal mass load is present (|Meiej[2003h . Next, we esti- 
mate the magnetization parameter of the jet plasma at the 
launching points, a, which reflects the effect of a rotating 
magnetic field on accelerating the jet plasma by measuring 
the Poynting flux in terms of particle flux (e.g. . iMichelll 19691 : 
ICamenzindii 19861 : iFendt fc GreineJl200lh . The initial mag- 
netization parameter of the jets (denoted by an index '0') is 
given by 



47rMjctsc3 



(36) 



where 'i! = j B dS and Qu is here taken as the angular 
velocity of the magnetic fleld lines frozen in the disc. For il- 
lustration, we evaluate ao for three values of the mass accre- 
tion rate (m = 1, 0.1, 0.01 and 0.001) in the case of a BH with 
the spin parameter of a, — 0.99 (see Fig. [S]). The magnetiza- 
tion parameter increases with decreasing mass accretion rate 
((To ~ rn^^), as well as with decreasing radius. When ao > 1, 
the Poynting flux dominates in the jets and the energy can 
be transferred from the magnetic field to the particles. As 
a result, the jets can be accelerated on the expense of the 



stored energy in the magnetic field as the Poynting flux is 
converted i nto kinetic energy flux by magnetic forces (e.g., 
lFendt|[2004h . The magnetic force can be split in two compo- 
nents: the magnetic pressure force (~ VB^ ), which points 
in positive outward direction, and the magnetic tension force 
(^ ), which points in negative inward direction. {B^ de- 
notes the strength of the toroidal component of the magnetic 
field.) In the case of a ballistic jet, which expands with a con- 
stant speed, the two forces cancel each other and the toroidal 
component of the magnetic fleld decreases as B^ ~ due 
to the magnetic flux conservation, {z denotes the distance 
along the jet.) To accelerate the flow, the magnetic pres- 
sure gradient must prevail over the magnetic tension force. 
This can be possible, in principle, due to the decrease in the 
strength of the toroidal magnetic field as the jet propagates 
away from the source, but one has to solve the full problem 
to verify whether the magnetic pressure gradient dominates 
over the magnetic tension force. The magnetic acceleration 
process is limited by the free energy available in the mag- 
netic field, and saturation must occur at some point. When 
the kinetic energy flux becomes dominant in the jet, strong 
shocks can occur. These shocks can further accelerate the 
jet. However, if the magnetic field becomes highly twisted, 
the magnetic field itself will not be able to explain the ac- 
celeration of the jet, as the magnetic reconnection or other 
instability in the jet can lead to magnetic energy dissipa- 
tion and shock formation. (Strong toroidal magnetic fields 
are subject to the kink instability, which excites large-scale 
helical motions that can distort or even disrupt the jet [e.g., 
iMizuno et al.l |200^ and references therein]. However, the 
growth rate of the kink mode may be reduced, for instance, 
by increasing the magnetic pinch or by including a gradual 
shear, an external wind or relativistic bulk motion.) How far 
the magnetic acceleration of the jets can occur depends from 
one case to another, and the ambient medium into which 
jets propagate can play a significant role. On pars ec scales, 
relati vistic shocks are expected to be prominent (jLobanovl 
|2007t ). A first stationary, strong shock can be produced in 
the approximate range (3 - 6) X 10^ rg (jMarkoff et al.ll200ll : 
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Figure 7. Lorentz factor of the jets as a function of m (Eq. 
I37|l . The jets have a bulk Lorentz factor 7 > 2 when the mass 
accretion rate m < 10~^'*, which corresponds to the spin-down 
power regime. In the case of the accretion power regime, m > 
10~^ *, the jets are mildly- and sub-rclativistic, 7 < 2. 



iMarscher et al.l[2008l ) , whereas moving shocks can occur be- 
tween 20 and 200 Tg. For very large values of the magne- 
tization parameter, the MHD approximation breaks down 
(jMizuno et al.ll2009l 'l. For (To = 1 the Poynting flux and the 
particle kinetic energy flux are in equipartition, whereas for 
ctq < 1 the jets are kinetic energy flux-dominated outflows. 
For given mass accretion rate and BH spin parameter, the 
value of (To is not much less than one, being only within one 
order or magnitude smaller than one. Therefore, it is possi- 
ble that the jets will be stable and propagate initially with 
a constant speed. 

For given mass accretion rate and radius, ao increases 
with the BH spin parameter; i.e., the Poynting flux in the 
jets increases as the BH rotates faster. For rh = 1 (Fig. [B]), 
the maximum value of ao is 0.993, which is obtained at the 
innermost stable orbit for a BH spin parameter of 0.9982. 
(However, ao increases for higher spins but we limit here the 
value of the BH spin parameter to that of Thome's model.) 
For rh = 10~^'*, which delimits the accretion power regime 
from the spin-down power regime, the initial magnetization 
parameter is larger than ~ 5 for any a, > 0.95. 

One can consider a relation between the magnetization 
parameter (Eq. I36|l and the initial Lorentz factor of the jets 
of the form 70 = af/'' , where the value of the power low index 
depends on the magnetic fleld conflguration . For a radia l 
outflow with negligible gas pressure, q = 3 (|Michellll969l ). 
For a collimated MHD je t, the value of q is also 3 if the flux 
distribution is the same (|Fendt fc Greinerll200ll '). With the 
power low above, the jets can present a radial distribution 
of the initial Lorentz factor, which increases with decreasing 
radius. Note that 70 differs from the bulk Lorentz factor 
of the jets (7) that is deflned below (Eq. I37|l . where it is 
assumed that the Poynting flux in the jets has been fully 
converted into kinetic energy flux. 

The bulk Lorentz factor of the jets 7, defined by 

-Pjcts = 7MjctsC^ = 7gjctsmA/EddC^, (37) 
[which follows, e . g., fro m lFalcke fc BiermannI (| 19951 ): see also 



IVila fc Romerol (|2010D ] is drawn in Fig. [7] as a function of 
the mass accretion rate. The jets have a relativistic speed of 
0.9 — 0.995 c (or 7 = 2—10, which is the typical bulk Lorentz 
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Figure 8. Rate of the disc angular momentum removed by the 
jets (Eq. I41|l as a function of the BH spin parameter a, for four 
values of the mass accretion rate. The mass accretion rates from 
top to bottom are 1, 0.5, 0.1 and 0.01. The bottom curve repre- 
sents the BH spin-down transferred to the disc inside the ergo- 
sphere by magnetic connection. Note that there is a slight differ- 
ence between the two bottom curves (red and turquoise curves). 



factor for an AGN jet) when the mass accretion rate rh £ 
-j^g-LSj. J these jets correspond to the spin-down 
power regime. In the case of the accretion power regime, 
the jets are mildly- and sub-relativistic (7 < 2). There is 
no significant variation of 7 with the BH spin parameter 
(a, ^ 0.95) for a given mass accretion rate. 



6 RATE OF THE DISC ANGULAR 

MOMENTUM REMOVED BY THE JETS 

Now, we define the rate of the disc angular momentum re- 
moved by the jets as 



2TvJL^rdr. 



(38) 



Using the angular momentum conservation law (Eg. llOf) . the 
rate of the disc angular momentum removed by the jets can 
be written as 

Jjcts = (1 - gjots) Mdc (lI - lL) +4tv [ rHdr, (39) 



where the dimensionless specific angular momentum of the 
gas particle orbiting in the accretion disc is 



, r, —3/2 , 2 —2 

_i/2 1 - 2a*r» +a«r* 



(1 - 3r^^ + 2a*rZ 



3/2 



1/2 ■ 



(40) 



Combining Eq. ^ with dHJ - (ETJ, (EHJ and (gTI, 
the rate of the disc angular momentum removed by the jets 
becomes 



Jjcts = rhMl^^cr^e 



Edd 



^(i-...)(4.-4.)(i^) 



_B_ 



V 109 Me 



ri-"Ry' (Oh. - t^D J dr. 



(41) 
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Figure 9. Efficiency of jet launching (Eq. 1 431 1 as a function of 
the BH spin parameter a» for different mass accretion rates. The 
mass accretion rates from the top to bottom are 0.0001, 0.001, 
0.01, 0.01 and 1. For very low mass accretion rates, rh < 0.001, the 
BH spin-down becomes a very efficient mechanism of launching 
the jets via the accretion disc. 



where C* and Rt are defined by Eq. L^^ and are 
the specific angular momentum of the gas particle (Eq. I40|) 
evaluated at the stationary limit surface and at the inner- 
most stable orbit, respectively. We can also write the disc 
angular momentum removed by the jets as the sum of two 
components, the accretion and the rotation parts. 



T Trace . tK 



(42) 



Figure |8] shows the rate of the disc angular momentum 
removed by the jets as a function of the spin parameter of 
the BH, given four values of the mass accretion rate (m — 1, 
0.5, 0.1 and 0.01), as well as the BH spin-down Jj°^^ (bottom 
curve) . To know how this angular momentum is transported 
by the jets, further models must be employed. 



7 EFFICIENCY OF JET LAUNCHING 

We define the efficiency of jet launching as the ratio of the 
launching power of the jets to the total power that comes 
from the gravitational energy of the accretion disc and from 
the rotational energy of the BH. Thus, 



V - 



mMEddc2 + P 



(43) 



jets 



In Fig. [9] we plot the efficiency of jet launching for the 
range of the mass accretion rate rh £ [0.0001, 1]. For very 
low mass accretion rates, rh < 0.001, the efficiency of jet 
launching reaches values close to unity, in which case the 
spin-down of the BH becomes a very efficient mechanism 
to launch the jets via the accretion disc. For the spin-down 
power regime, the efficiency of jet launching is higher than 
the maximal efficiency of converting the gravitational energy 
of the accretion disc into radiation (e = 0.42), as a result of 
transferring the BH rotational energy to the accretion disc 
via BH-disc magnetic connection. 



8 SPIN EVOLUTION OF THE BLACK HOLE 

Theoretically, a Kerr BH can be spun up to a state with 
a spin parameter whose maximum value is a* = 1. As the 
spin evolves, a Kerr B H can achieve a stationary state. A 
theorem established bv lHawkind j 19721 ) states that a BH is 
in a stationary state if and only if the BH is either static 
or axisymmetric. Suppose we have a Kerr BH. Perturbing 
fields can, however, defied the spin orientation away from 
the symmetry axis. In this case, the BH must either spin 
down until a static (Schwarzschild) BH is reached or evolve 
in such a way that it aligns its spin with the perturbative 
field orientation. 

Next, we study the BH spin evolution and seek the max- 
imum spin parameter that corresponds to a stationary state 
of the BH, when both the BH-dis c magnet i c con nection and 
the jet formation are considered. iThornd (|l974l ) calculated 
the infiuence of photon capture on the spin evolution of the 
BH and found a limiting state of a^jim — 0.9982. This limit 
does not apply to the model proposed here since the disc 
inside the BH ergosphere is not radiant, as in the case of 
Thome's model. Instead, it drives the jets. We consider this 
limit only to determine the maximum value of the BH mag- 
netic field, given at the time when the BH accretes at near 
the Eddington limi t. 

iBardeenI (|l970l ) showed that the mass and angular mo- 
mentum of the BH can be changed by the specific energy and 
angular momentum of the particles carried into the BH. The 
BH mass (and the angular momentum) variation equals the 
value of the particle specific energy (and angular momen- 
tum) at the innermost stable orbit multiplied by the rest 
mass accreted (dMo) if no other stress energy is allowed to 
cross the horizon. That is. 



dM = Elns dMo and dj = J^. dMo 



(44) 



2a.. 



(45) 



where i?^a Skud are the specific energy and angular mo- 
mentum of the particles evaluated at the innermost stable 
orbit. Using Eq. (|44|) . one can obtain the differential equa- 
tion that describes the spin evolution of the BH due to mat- 
ter accretion: 

da, \ c / dJ \ 

dlnM J "gmIZmJ 

\ / matter \ / 

Now, we consider the magnetic extraction of the BH 
rotational energy through the BH-disc magnetic connection. 
The spin evolution law (Eq. I45p will be changed due to the 
counter-acting torque exerted on the BH by the magnetic 
field that connects the BH to the disc inside the ergosphere. 
The energy and angular momentum lost (or gained, depend- 
ing on the angular velocities of the BH and disc, cf. Eg. I14p 
by the BH th rough the BH-disc magnetic connection are 
(e.g..[Lil l2002al l: 

(^) = 2Phd and (^) = 2Thd, (46) 

\ / HD \ / HD 

where Phd ~ ^hTud, and the factor '2' comes from the fact 
that the accretion disc has two surfaces. Adding the effects 
of the BH spin-up by accretion (Eq. I44p and the BH spin- 
down by magnetic connection (Eq. I46[) . the equations for 
evolution of the BH mass and the BH angular momentum 
become: 



fdM\ 



= (1 - fflets 



) Mj,cEI, 



+ C 



(dM\ 



(47) 
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Figure 10. Spin evolution of a Kerr BH. The purple line repre- 
sents the driving torque by which the matter spins up the BH (Eq. 
I45I I. The blue and green lines represent the counter-acting torque 
on the BH due to transfer of rotational energy from the BH to the 
disc (Eq. I50I I for m = 0.1 and m = 0.01, respectively. For these 
two cases of the mass accretion rate, the stationary state of the 
BH corresponds to a» ~ 0.9944 and a» ~ 0.9762, respectively. 



(48) 



Using these two equations, we can express the BH spin evo- 
lution as 



da* 



d In M 



C (1 - gjcts)A/DiL + 



dj_\ 

dt / hp 



2a., 



(49) 

when both the BH spin-up by accreting matter and the BH 
spin-down due to the angular momentum transferred from 
the BH to the disc inside the ergosphere are considered. 
From Eqs. (|45|l and (|49p . the spin-down of the BH by means 
of BH-disc magnetic connection is described by 



d ]nM 



da* 
d \nM 



dat 
d InM 



(50) 



Figure [To] shows the spin evolution of a Kerr BH for two 
values of the mass accretion rate m = 0.1 and m = 0.01, re- 
spectively. The purple line represents the driving torque by 
w hich the matte r spins up the BH [Eq. 1451 see also fig. 6 
m iThornd (|l974h ]. and the blue and green lines represent 
the counter-acting torque on the BH due to transfer of rota- 
tional energy from the BH to the disc (Eq. I50p for m = 0.1 
and m = 0.01, respectively. The crossing point of the plots 
corresponds to the spin parameter for which the BH is in 
a stationary state. For m = 0.1, if the BH initially rotates 
with a* = 0.9982, the BH may spin down to a stationary 
state with a maximum spin parameter of a, = 0.9944. For 
rh = 0.01, the maximum spin parameter is a* = 0.9762, 
whereas for m = 0.001, the maximum spin parameter is 
a* = 0.9525. Thus, as the mass accretion rate decreases, the 
maximum spin parameter corresponding to a stationary BH 
decreases as well. On the other hand, as the mass accretion 
rate decreases, the magnetic torque reaches values close to 
unity, so it is greater than 0.43, which is the maximum value 
of the matter torque. This implies deviations from pure Ke- 
plerian orbits, and so the possibility to drive away the excess 
angular momentum of the disc in the form of jets when the 
BH spin-down power is considered. A further analysis of 



the spin evolution (which is not explicitly shown in Fig. I10|l 
suggests that a BH needs a mass accretion rate of at least 
m ~ 0.001 for its spin to stay high (a* ^ 0.95). For lower 
mass accretion rates (m < 0.001), the BH may spin down 
continuously until the BH reaches a static state. It can spin- 
up again to a* ^ 0.95 if a large amount of matter is provided 
by accretion (or by merging, which is not discussed here). In 
this case, the amount of accreted mass shoul d be a factor of 
about 1.84 from the initial mass of the BH (|Thornelll97i ). 



9 MAXIMUM LIFETIME OF THE AGN FROM 
THE BH SPIN-DOWN POWER 

In this section, we calculate the time-scale needed for a Kerr 
BH to spin down from a* ~ 1 to 0.95, which can then be 
related to the maximum lifetime of the AGN, provided that 
the BH was spun up to nearly its maximum spin during a 
phase when the AGN was active. The AGN can be active 
as long as the gravitational energy of the accretion disc is 
converted into observed radiation energy. Such an AGN can 
have a longer lifetime through the additional use of its BH 
spin-down power, despite having a very-low mass- accretion 
rate ( m < 0.1). Following the well-known work by ISalpeterl 
(|l964h . the time needed to fuel the AGN to a bolometric 
luminosity Lboi ~ 10*^ erg s"'^ can be ~ 10^ yr for a typical 



radiative efficiency of e = 0.1. Moreover, the lifetime of high 
accreting AGN (and quasars) was constraine d by recent ob- 
servations to the range ~ 10^ — 10^ yr (e.g.. [Porciani et al.l 
|2004 iHopkins fc Hernauist|[2009l ) . 

Next, we estimate the maximum lifetime of the AGN. 
Differentiating the BH angular momentum J = Mca = 
{GM'^ /c) a, with respect to time t, the BH time evolution 
is specified by 

Integrating this equation, the time interval over which the 
BH spin evolves between two given values of a, is 



da, 
~dt 



t 



GAP 



a, fdM\ 
M \ dt ) 



da,, (52) 



where (dj/dt) can be obtained from Eqs. (|48p and (|50p . and 
(dM/dt) from Eqs. (|47)) and ((50)) . With the above equation, 
we can estimate the lifetime of the AGN. The time interval 
is not dependent on the BH mass. 

In Fig. 111! we plot the time evolution of the AGN as 
a function of the mass accretion rate (Eq. [52}, when the 
BH spin parameter decreases from a, '-^ 1 to 0.95. For a 
mass accretion rate close to the Eddington limit, the life- 
time of the AGN is about 3 x 10^ yr. The lifetime curve 
moves toward lower mass accretion rates for another ~ 10* 
yr, when the AGN uses its BH spin-down power to launch 
the jets. Therefore, the total lifetime of the AGN can be 
much longer than the Hubble time {tn ~ 10^"'^* yr). The 
maximum lifetime of the AGN is, however, dependent on 
the mass accretion rate. The maximum lifetime of the AGN 
from the BH spin-down power is, for instance, ~ 2.8 x 10** yr, 
~ 3.9 X 10* yr and ~ 4 x 10* yr for m = 10"^ m = 10~'^ and 
m < 10"*, respectively. In the latter case, the BH may not 
attain a stationary state and spins down until a static BH 
is reached. The lifetime of the AGN scales with the strength 
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Figure 11. Lifetime of the AGN from the BH spin-down power as 
a function of the mass accretion rate. The hfetime of the AGN is 
~ 3 X 10^ yr when using the accretion power. The BH spin-down 
power adds to the hfetime of the AGN, for instance, 2.8 X 10* yr, 
3.9 X 10* yr and ~ 4 X 10* yr when m = 10~^, ih = 10^^ and 
rh < 10~*, respectively (while the BH spin decreases from a, ~ 1 
to 0.95). 



of the BH magnetic field relative to its maximum value as 
t ~ (Bh/Bh'''')'^. Therefore, if the Bh is a factor of k lower 
than Bh'^") the maximum lifetime of the AGN will be a fac- 
tor of A:^ larger. For instance, when k = 7 and rh = 10~^, 
one obtains the exact Hubble time. 

Now, we compare our results to the lifetime of an AGN 
powered by the Blandford-Znajek mechanism. The total en- 
ergy that can be extrac ted by the Blandford-Znajek mech- 
anism is (e.g..lLi| [2000al '): 

Ebz ^ 0.09Mc^ ~ 1.6 X 10"' erg (j^j^^ ■ (53) 

The maximum rate of energy ext raction by the Blandford- 
Znajek mechanism is [Eq. 4.50 of lThorne et all (|l986l )]: 

for a BH with a* ~ 1. Therefore, the lifetime of an AGN 
powered by the Blandford-Znajek mechanism is: 

tf = ^ ifR;)"(T5^)"' <^^' 

which is not dependent of the BH mass, as Bh scales with 
{M/10^ Mq)~^^'^ . Except for the exact time-scale, our re- 
sult (Eq. I52p scales with the magnetic field and with the 
BH mass exactly the same way as for the Blandford-Znajek 
mechanism. 

Moving back to our results, if the mass accretion rate 
changes over the whole life of the AGN from m ^ 10~^ * to 
m < 10~*, the maximum lifetime of the AGN can be even 
longer than that of an AGN powered by the Blandford- 
Znajek mechanism. In summary, the maximum lifetime of 
the AGN can be much longer than ~ 10*^ yr when using 
the BH spin-down power. The lifetime is dependent on the 
mass accretion rate, as well as on the factor (Bh/Bh^")- 
We mention that the results presented here refer only to 
rapidly-spinning BHs (a, ^ 0.95). 



10 SUMMARY AND CONCLUSIONS 

Starting from the general-relativistic conservation laws for 
matter in a thin accretion disc, we included both the BH- 
disc magnetic connection and the jet formation. The jets are 
launched from the disc inside the BH ergosphere, where the 
rotational effects of the space time become much stronger. In 
the BH ergosphere, the frame-dragging effect ensures that 
the magnetic field lines frozen in the disc co-rotate every- 
where with the BH, reducing the effect of magnetic recon- 
nection and keeping the magnetic field configuration glob- 
ally the same. Furthermore, the jets can extract mass, en- 
ergy and angular momentum from the inside the ergosphere 
disc. For this situation, we derived the mass flow rate into 
the jets, the launching power of the jets, the angular mo- 
mentum removed by the jets, the efficiency of launching the 
jets, the maximum spin parameter attained by a stationary 
BH and the maximum lifetime of an AGN. 

• We found that the mass flow rate into the jets is depen- 
dent on the BH spin parameter, where the mass outfiow into 
the jets can be associated with the rotation of the space-time 
itself. For near maximal spin parameter a* ~ 1, the mass 
flow rate into the jets is about 98 per cent of the available 
mass flow in the disc inside the BH ergosphere, whereas for 
a* = 0.95 this is only about 8 percent. This means that in 
the case of maximal spin, the BH almost stops being fed by 
accreting matter. As a possible alternative, the jets may have 
no matter right at the beginning (i.e., they are Poynting flux 
jets) and get it only very quickly from drifts just above the 
disc or the surrounding wind (i.e., indirectly from the disc). 
This may only be relevant for extremely low accretion rates. 

• In this work, we considered the case of rapidly-spinning 
BHs with a spin parameter of a* ^ 0.95 and a mass of 
IO'-'A/q, and we assumed that the power of the disc inside the 
BH ergosphere is used to drive the jet. We found that at low 
mass accretion rates, the jet power can be supplied by the 
BH rotational energy via the disc inside the ergosphere. The 
switch from an accretion power regime to a spin-down power 
regime corresponds to a mass accretion rate of m ~ 10~^ *. 
In the case of the spin-down power regime (m < 10"^ *), 
the jet power is ^ 10*'' erg s~^, which is only 10"' of the 
Eddington luminosity of a 10^ solar mass BH. This is com- 
parable to the maximum rate of energy extraction by the 
Blandford-Znajek mechanism, which is ~ 10^^ erg s~^ for a 
BH mass of 10^ Mq and a, ~ 1. This implies that, in prin- 
ciple, both the Blandford-Znajek mechanism and launching 
jets from the disc inside the ergosphere via the BH-disc mag- 
netic connection can operate. We intend to study the driving 
of the jets when the BH is threaded by a combination of open 
and closed magnetic held lines in future work. 

• The jets can have a relativistic speed, 0.9 — 0.995 c (or 
7 = 2 — 10, which is the typical bulk Lorentz factor for an 
AGN jet), when the mass accretion rate m € [10~''^, 10~^ *]. 
In the case of the accretion power regime, the jets are mildly- 
and sub-relativistic. However, after launching the jets can 
be accelerated through magnetic processes. The jets remove 
the angular momentum of the disc inside the BH ergosphere 
at a rate which is dependent on the BH spin parameter. 
To know how this angular momentum is transported by the 
jets, further models have to be employed. The efficiency of 
jet launching is higher at low mass accretion rates, reaching 
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values close to unity for m ~ 10^'*. In this case, the BH 
spin-down power is efficiently used to launch the jets. 

• Considering the balance between the BH spin-up by ac- 
creting matter and the BH spin-down due to the magnetic 
counter-acting torque on the BH, we determined the maxi- 
mum spin parameter which corresponds to a BH stationary 
state. The maximum spin value shifts towards a» — 0.95 as 
the mass accretion rate decreases. For instance, the maxi- 
mum spin parameter corresponding to m = 0.1, rh = 0.01 
and m = 0.001 is a, = 0.9944, a. = 0.9762 and a. = 0.9525, 
respectively. At lower mass accretion rates (m < 0.001), 
the BH may undergo a spin-down process towards a static 
BH. The BH never reaches a stationary state unless a large 
amount of matter is provided (perhaps by star capture or 
by merging) to spin up the BH again to a, ^ 0.95. 

• We showed that an AGN can have a much longer life- 
time than ~ lO'^ yr when using the BH spin-down power, 
and the maximum lifetime is dependent on the mass ac- 
cretion rate, as well as on the factor (Bh/Bh^"). After an 
accretion-dominated phase of about 3 x 10^ yr, the AGN 
can live off of the BH spin-down power for another 10* yr. 
The BH spin-down power adds to the lifetime of the AGN, 
for instance, ~ 2.8 x 10® yr, ~ 3.9 x 10* yr and ~ 4 x 10* yr 
for rh — 10~^, m — 10^^ and m < 10~*, respectively. More- 
over, if the Bn is a factor of k lower that B^^^, then the 
lifetime of the AGN will be a factor of fc^ larger. For k = 7 
and m = 10~^, one obtains the exact Hubble time. Another 
possibility is that the mass accretion rate changes over the 
whole life of the AGN from rh ~ 10"^ * to m < 10"". In 
this case, the maximum lifetime of the AGN can be even 
longer than that of an AGN powered by the Blandford- 
Znajek mechanism, which is ~ 5 x lO" yr for a BH with a, 
close to the maximal value. However, it will be difficult to 
predict a maximum lifetime of the AGN for this possibility, 
since there is no mechanism, to date, to control the change 
of the mass accretion rate over long intervals of time. 

In the limit of the BH spin-down regime, the model pro- 
posed here can be regarded as a variant of the Blandford- 
Znajek mechanism. It is a variant of the Blandford-Znajek 
mechanism if an insignificant fraction of the rotational en- 
ergy is transported away by jets through open magnetic 
field lines that thread the BH. If we were to consider the 
open magnetic field lines in addition to the closed magnetic 
field lines that thread the BH, our calculations would have 
been completed such that they would include the Blandford- 
Znajek mechanism. 

The results presented in this paper are dependent on 
our assumptions that a BH-disc magnetic connection exists. 
Closed magnetic field lines in the BH ergosphere may be pro- 
duced by a current ring in the vicinity of the BH. Models 
for the magnetic connection where a poloidal magnetic field 
is generated by a single electric current flowing in the BH 
equatorial plane or at the inner ed ge of t he a ccretion disc 
were p roposed, for instance, byQ (|2002 j ) and lWang et al.l 
(I2OO3). The key parameters of the proposed model are, how- 
ever, the BH mass, the BH spin and the mass accretion rate. 

One indirect way to test whether this mechanism oper- 
ates in reality is to study the relation between the observed 
radio flux density from AGN (e.g., from flat-spectrum core 
sources) and their mass accretion rates in order to fit the 
model prediction with respect to the relation between the 



power of the jet and the mass accretion rate (see Fig. |4]). 
That relation shows that the power of the jet does not de- 
pend linearly on the mass accretion rate all the way down to 
very low accretion rates, so that there can be sources with 
relatively strong jet power but low mass accretion rate. In 
this case the jet power is mainly dependent on the BH pa- 
rameters, such as the mass and the spin of the BH. The 
first step is to find a relation between the power of the jets 
and the observed radio flux density, -Fobs- For a conical jet 
(where the magnetic field along the jet scales with the dis- 
tance, B ~ 2~^, and the electron number density in the 
jet scales as C' ~ 2~^) from a flat-spectrum core source, 
we found the dependence: Pjots ~ F^^bC^D^^/^M"''/", where 
D is the distance to the AGN and M is the BH mass. A 
full description of its derivation is the subject to a paper 
in preparation. The second step is to produce a complete 
sample of AGN with known jet parameters, as the Doppler 
factor, and whose mass accretion rate can be constrained by 
observational data, and then to fit the model prediction. 

Furthermore, if AGN were in the Poynting flux limit 
disregarding the mass accretion rate, then one might expect 
that the power in the jet from a large sample runs into a 
lower limit, which is given by the minimum Poynting flux. 

The model presented here can also be extended to mi- 
croquasars assuming that physical quantities (e.g., BH mag- 
netic field, jet power, etc.) scale with the BH mass. 

Although numerical simulations of jet formation from 
the ergosphere of a rapidly-spinning BH with closed mag- 
netic field lines that connect the BH to the disc inside the 
ergosphere has not been performed yet, this can be one of 
the challenges to be faced by numerical relativists. 
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APPENDIX A: DERIVATION OF EQUATION 4 

Here, we present the derivation of the angular momentum 
conservation law for the matter that flows in the accre- 
tion disc (Eq. I10|l . This derivation is based on the general- 
relativistic angular-momentum-conservation law that de- 
scribes the structure o f a geometrically thin accretion disc 
Page fc Thornelll974 l. Based on this conservation law, iLiI 



2002al ) derived the conservatio n law that includes t h e BH - 
disc magnetic connection, and iDonea fc BiermannI (|l996l ) 
derived the conservation law that includes the jet. In a 
slightly different manner, our result is obtained when both 
the BH-disc magnetic connection and the jet formation are 
considered. 

The angular-momentum conservation law of the matter 
in a standard, t hin a ccretion disc is given by eq. (23) of 
IPage fc Thornd i I974I ). For details on this procedure and 
the meaning of the notations used, t he reader is referred 
to the paper bv lPage fc Thornd (|l974 ). Below, we consider 
the first and second integrals of their equation, keeping the 
non-vanishing terms: 



r + Ar 



(Al) 



(27rAt) (< po > U4,u+ < >) \f\g\dz \ 

= {(27rAt)(ELV + W^;)\^}^^^'' 

= At Ar [- (Md - Mjets + 27rrlV;] , 

where we made use of the rest-mass conservation law that 
includes the mass fiow into the jets , Md — Afjots = — 27rrEu'', 
and 

< / (27rAf) (< po > u^u'' < q~ >) \/\g\dr > 

= 2(27r Af)(Ew^ +F)L'^y^\Ar = 2{2Tir AtAr)JL\ 

(A2) 

where J = Eu' -I- F denotes the total flux of energy (of 
particle and electromagnetic origin) transported by jets. 

Adding the two integrals, eq. (23) of Page fc Thorne 
becomes: 

[(A/d - Afjets)L^ - 27rrVy;] = 4ttJL\ (A3) 

Now, using the definition of the magnetic torque pro- 
duced by the BH on both surfaces of the accretion disc which 
is given by Li (see Eqs. [iT] and 114(1 and including c, Eq. IA3I 
becomes: 

^ [(1 - gjets) A/dcL^I + 4TTrH = 4nrJL'', (A4) 

when both the BH-disc magnetic connection and the jet for- 
mation are considered. 



